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Next-to-leading order QCD corrections to e + e — > 3 jets 
with massive quarks 

A. Brandenburg*, W. Bernreuther, and P. Uwer 

Institut fur Theoretische Physik, RWTH Aachen, 52056 Aachen, Germany 

We discuss a calculation of the next-to-leading order QCD corrections to the process e + e~ — > 3 jets with 
massive quarks, and show numerical results for the three jet fraction and the differential two jet rate. 



1. INTRODUCTION 

In view of the large number of jet events at 
the Z resonance collected both at LEP and SLC, 
it is desirable for precision tests of the stan- 
dard model to use in the theoretical predictions 
next-to-leading order (NLO) partonic matrix ele- 
ments that include the full quark mass depen- 
dence. This is particularly important for b quark 
enriched samples, which can be obtained with 
high purity using vertex detectors. 

The NLO QCD corrections to e + e~ — > 3 jets 
are well-known for massless quarks The 
three, four, and five jet rates involving massive 
quarks have been computed to leading order in a s 
already some time ago ||,|| . Recently, results for 
the NLO corrections to e + e~ — > 3 jets including 
mass effects have been reported |7|-|TT[. Knowing 
the NLO matrix elements including the mass de- 
pendence, one can try to extract the mass of the 
b quark from three jet rates involving b quarks at 
the Z peak. This was suggested in |Q, elabo- 
rated in |Q,^,^3|, and experimentally pursued by 



the DELPHI collaboration 14 13]. Further appli 



cations include precision tests of the asymptotic 
freedom property of QCD by means of three jet 
rates and event shape variables measured at var- 
ious center-of-mass energies, also far below the Z 
resonance |l(|. For theoretical predictions con- 
cerning the production of top quark pairs at a fu- 
ture e + e~ collider, the inclusion of the full mass 
dependence is of course mandatory. 

We have computed the complete differential 
distributions for e + e~ annihilation into three and 



four partons via a virtual photon or Z boson at 
order a?, including the full quark mass depen- 
dence • This allows for order a 2 predictions 
of oriented three jet events (for the massless case, 
see [0), and of any quantity that gets contribu- 
tions only from three and four jet configurations. 
In this talk we give an outline of our calculation 
(for details, see [|ll|), after which we present some 
numerical results. 

2. OUTLINE OF THE CALCULATION 

The calculation of an arbitrary quantity domi- 
nated by three jet configurations and involving a 
massive quark-antiquark pair to order a 2 consists 
of two parts: First, the computation of the ampli- 



tude of the partonic reaction e 



j*,Z* 
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QQg at leading and next-to-leading order in the 
QCD coupling. Here Q denotes a massive quark 
and g a gluon. Second, the leading order ma- 
trix elements of the four-parton production pro- 
cesses e+e" -> Z*,7* -> ggQQ,QQqq,QQQQ 
are needed. Here q denote light quarks which are 
taken to be massless. 

The infrared (IR) and ultraviolet (UV) singu- 
larities, which are encountered in the computa- 
tion of the one-loop integrals, are treated within 
the framework of dimensional regularization in 
d = 4 — 2e space-time dimensions. We remove 
the UV singularities by the standard MS renor- 
malization. We have converted from the outset 
the on-shell mass of the heavy quark Q into the 
corresponding running MS mass. It is known that 
far from threshold one thereby absorbs some large 
logarithms into the running mass. 

After renormalization, the virtual corrections 



2 



to the differential cross section for e + e~ — > QQg 
still contain IR singularities. These have to be 
cancelled by the singularities that are obtained 
upon phase space integration of the squared tree 
amplitudes for the production of four partons. 
Different methods to perform this cancellation 
have been developed (see [^8|-p0| and references 
therein). We use the so-called phase space slic- 
ing method elaborated in jl8| , which wc modified 
to account for masses pl |. The basic idea is to 
"slice" the phase space of the four parton final 
state by introducing an unphysical parton resolu- 
tion parameter s m in, which is much smaller than 
all relevant physical scales. The parameter s m in 
splits the phase space into a region where all four 
partons are "resolved" and a region where at least 
one parton remains unresolved. In the unresolved 
region soft and collinear divergences reside, which 
have to be isolated explicitly to cancel the singu- 
larities of the virtual corrections. This is con- 
siderably simplified due to collinear and soft fac- 
torizations of the matrix elements which hold in 



the limit s r 



0. (In the presence of massive 



quarks, the structure of collinear and soft poles 
is completely different as compared to the mass- 
less case.) After having cancelled these IR poles 
against the IR poles of the one-loop integrals en- 
tering the virtual corrections, one is left with a 
completely regular differential three-parton cross 
section which depends on s m in ■ The contribution 
to a three jet quantity of the "resolved" part of 
the four-parton cross section is finite and may be 
evaluated in d = 4 dimensions, which is of great 
practical importance. It also depends on s m in 
and is most conveniently obtained by a numeri- 
cal integration. Since the parameter s m in is com- 
pletely arbitrary, the sum of all contributions to 
any observable must not depend on s m i n . Since 
the individual contributions depend logarithmi- 
cally on s m in, it is a nontrivial test of the cal- 
culation to demonstrate that observables become 
independent of s m in for small values of this pa- 
rameter. 

3. NUMERICAL RESULTS 

We will now discuss results for some three 
jet observables involving massive quarks. All 
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Figure 1. The three jet fraction /| at NLO as de- 
fined in the text at ^J~s = y, = mz as a function of 
y-min = Smin /{sy cut) for the Durham algorithm at a 
value of the jet resolution parameter y cut = 0.03 with 

m^ifi = m z ) = 3 GeV and a s (fj, = m z ) = 0.118. 

quantities are calculated by expanding in a s to 
NLO accuracy. We consider here the JADE 
pl| and Durham |22l jet finding algorithms, al- 
though other schemes [p3[ can also be easily im- 
plemented. First we demonstrate the indepen- 
dence of physical quantities on the parameter 
Smin as s m in —> 0. We choose as an example 
the three jet fraction for b quarks, 



ft{y cut) 



(1) 



the numerator u\ is defined as the three 



In Ql| 

jet cross section for events in which at least two 
jets containing a b or b quark remain after the 
clustering procedure. This requirement ensures 
that the cross section stays finite also in the 
limit m& — > 0. The contribution of the process 
e + e~ — ► 7* — ► qqg* — » qqbb to the three jet 
cross section with one tagged b quark develops 
large logarithms ln(m^) - which find no counter- 
part in the virtual corrections against which they 
can cancel - when the bb pair is clustered into a 
single jet. In principle, there are three ways to 
handle this problem: One may impose suitable 
experimental requirements/cuts to suppress con- 
tributions from events with two light quark jets 
and one jet containing a bb pair (the definition for 
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Figure 2. The three jet fraction /f as defined in the 
text at y/s = mz as a function of y cu t for the Durham 
algorithm. The dashed line is the LO result. The 
NLO results are for ii = mz (solid line), fi = mz/2 
(dotted line), and /i = 2m z (dash-dotted line). 
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Figure 3. Dependence of the three jet fraction /| for 
the Durham algorithm on the renormalization scale 
at y/s = mz and y cu t = 0.2 x 1CT 3/10 « O.f for on- 



shell b quark masses m^ olc = 3 GeV (full curves) and 



e = 5 GeV (dashed curves) 



(Tg chosen by us is an example for this), or one can 
improve the fixed order calculation by absorbing 
the large logarithm into a fragmentation function 
for a gluon into a b quark. The third possibility is 
simply to keep the large ln(mf) term. This may 
however lead to an overestimate of the b quark 
mass effects in some observables. A detailed dis- 
cussion of this issue will be presented elsewhere 
p3 . Fig. [j] shows the three jet fraction /| in the 
Durham scheme at NLO as a function of y m in = 
s min /(sy cu t) at a fixed value of y cut = 0.03 and 
y/s = mz = 91.187 GeV. The corresponding plot 
for the JADE scheme has a similar shape. For the 
renormalization scale we take in this plot /i = y/s. 

As mentioned above, we use (/i) defined in 
the MS scheme at the scale //. The asymptotic 
freedom property of QCD predicts that this mass 
parameter decreases when being evaluated at a 
higher scale. (For low energy determinations of 
the b quark mass see e.g. 25|2^| and references 
therein.) With mf^ifi = m b ) = 4.36 GeV [|| 
and a s (mz) — 0.118 as an input and employ- 
ing the standard renormalization group evolution 
of the coupling and the quark masses, we use 
the value 



if s (fj, = m z ) = 3 GeV. One clearly 
/| reaches a plateau for small values 
The error in the numerical integration 



and y min 



becomes bigger as y m in — > 0. In order to keep 
this error as small as possible without introduc- 
ing a systematic error from using the soft and 
collinear approximations, we take in the follow- 
= 0.5 x 10~ 2 for the Durham algorithm 
10~ 2 for the JADE algorithm. At 
these values, the dominant s min -dependent indi- 
vidual contributions from three and four resolved 
partons are about a factor of 4 (Durham) and 
2.5 (JADE) larger than the sum. In Fig. || we 
plot /| as a function of y cu t, again at y/s = mz- 
The QCD corrections to the LO result are quite 
sizable as known also in the massless case. The 
renormalization scale dependence (where /i is var- 
ied between mz/2 and 2m z), which is also shown 
in Fig. H is modest in the whole y cu t range exhib- 
ited for the Durham algorithm. In Fig. || we take 
a closer look on the scale dependence of /|, now 
using the on-shell mass renormalization scheme. 
We vary the scale fj. between mz/16 and 2mz for 
a fixed value y cut = 0.2 x 10~ 3 / 10 « 0.1 and on- 



shell masses w^ ole 



3 GeV and m 



pole 



5 GeV. 



We see that the scale dependence of the LO result 
(which is solely due to the scale dependence of a s 
at this order) in the Durham algorithm amounts 
to about 100% in the /x interval shown. The in- 
clusion of the at corrections reduces the scale de- 
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Figure 4. The double ratio C as a function of the 
cm. energy at y cu t = 0.08 for the JADE algorithm. 
The full curve and the points show the NLO result, 
the dashed curve shows the LO result. In both cases a 



running mass m^ S 



(fj, = y/s) evolved from m b (fJ> = 
mz) = 3 GeV is used. The dash-dotted curve shows 
the LO result for a fixed mass m/, = 4.7 GeV. 



pcndcnce significantly; the NLO result for /| at 
/i = 2mz is about 10% smaller than the NLO 
result at [i = m^/16. In the case of the JADE 
algorithm, which we do not show here, the differ- 
ence between /| at /i = 2mz and at n = mz/16 
is reduced from about 100% at LO to about 30% 
at NLO. 

The effect of the b quark mass may be illus- 
trated by looking at the double ratio 



C{y cu t) 



(2) 



where the denominator is the three jet frac- 
tion when summing over all active quark flavors, 
which is given to a very good approximation by 
the massless NLO result f27|j . Similar double ra- 
tios have been studied in |7j-^], and |l^ , [l3|| . In 
Fig. H we plot C as a function of the cm. en- 
ergy at y cut = 0.08 for the JADE algorithm. The 
running of a s is taken into account in the curves, 
where we again use as an input a s (fi = mz) = 
0.118. For the dashed (LO) and full (NLO) curve 
we use a running mass (/x = y/s) evolved 
from m^ IS (/x = mz) = 3 GeV. For comparison 
we also show the LO result for a fixed value of 
the b quark mass m& = 4.7 GeV (dash-dotted 



Fi gure 5. The double ratio T> as a function of y cu t 
for y/s = n = mz- Full circles: NLO results for 
m^ /[s (^i = mz) = 3 GeV. For comparison, the squares 
(triangles) are the LO results for mj = 3 GeV (mj = 
5 GeV) . The horizontal bars show the size of the bins 
in y C ut- 

curve), which is the corresponding value of the 
pole mass. One clearly sees that the effect of the 
b quark mass gets larger for smaller cm. energies. 

Another interesting quantity to study mass ef- 
fects is the differential two jet rate Q defined 
as 



D 2 (y) 



h{v) - h{y- &y) 

Ay 



(3) 



where /a(y) is the two jet fraction at y — y cu t for 
a given jet algorithm. The advantage of D2 over 
the three jet fraction fa lies in the fact that the 
statistical errors in bins of D2(y) are independent 
from each other since each event enters the dis- 
tribution only once. To order oq, Z?2 (j/) can be 
calculated from the three- and four jet fractions 
using the identity 



1 = h + h + h + 0(a 3 s ). 
We define 



2%) 



D\{y) 



nincl 



(y) 



(4) 



(5) 



where we - as in the case of the quantity C - use 
the massless NLO result to evaluate the denomi- 
nator. We plot our result for T>(y) in the Durham 



5 



scheme in Fig. ^, again for y/s = fi = mz- The 
effects of the b quark mass are of the order of 5% 
or larger at small values of y cu t- 

4. CONCLUSIONS 

We have presented NLO results for a number 
of three jet observables involving massive quarks. 
It will be interesting to see how our predictions 
compare to detailed experimental analyses. 
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DISCUSSIONS 
A. P. Contogouris 

You stated three ways for treating large logarithms 
and that you followed one of them. Suppose you 
follow any of the other two. How much your 
higher corrections will change? Do you have any 
idea about their sensitivity, perhaps from another 
process ? 

A. Brandenburg 

// one adds the contribution from g* — > bb split- 
ting "naively", i.e. keeps the large logarithm in 
the three jet cross section with one tagged b quark 
as it is, the ratio C(y cu t) evaluated at i/s = [i = 



G 



mz is, for large values of y cu t, U P t° 3% larger 
as compared to the result for this ratio with two 
tagged b quark jets. The difference decreases for 
smaller values ofy cu t- 



